RPAE versus RPA for the Tomonaga model with quadratic energy dispersion 
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Recently the damping of the collective charge (and spin) modes of interacting fermions in one 
spatial dimension was studied. It results from the nonlinear correction to the energy dispersion in 
the vicinity of the Fermi points. To investigate the damping one has to replace the random phase 
approximation (RPA) bare bubble by a sum of more complicated diagrams. It is shown here that 
a better starting point than the bare RPA is to use the (conserving) linearized time dependent 
Hartree-Fock equations, i.e. to perform a random phase approximation (with) exchange (RPAE) 
calculation. It is shown that the RPAE equation can be solved analytically for the special form 
of the two-body interaction often used in the Luttinger liquid framework. While (bare) RPA and 
RPAE agree for the case of a strictly linear disperson there are qualitative differences for the case 
of the usual nonrelativistic quadratic dispersion. 



quires the solution of an integral equation£i£ 

For spin 1/2 fermions collective modes usually occur 
in the charge and the spin density response. They were 
studied by Dzyaloshinski and Larkin^ using standard 
many body techniques (avoiding bosonzation). They 
showed that for linear dispersion the exact results for the 
collective modes for the Tomonaga model are identical to 
the RPA result. Diagramatically this can be understood 
in terms of the "closed loop theorem" which holds if one 
assumes strictly linear energy dispersion around the two 
Fermi points^. 

Recently various activities started to study the charge 
and spin density response when deviations from a strictly 
linear dispersion for one-dimensional fermions is taken 
into accoun t 10 ! 11 ' 12 ! 13 . Then the approximation using the 
"bare bubble" is no longer exact. 

It is one purpose of this paper to point out that the 
bare RPA is not even a good starting point for a pertur- 
bative calculation as e.g. it has the one-particle one-hole 
(lp-lh) continuum in the wrong place. It is much better 



I. INTRODUCTION 

In a seminal paper Tomonaga 1 presented the exact so- 
lution for the long wavelength density response of a sys- 
tem of interacting fermions in one spatial dimension. In 
order to simplify the problem Tomonaga studied the high 
density limit where the range of the interaction is much 
larger than the interparticle distance. Then the Fourier 
transform v(k) of the two-body interaction is nonzero 
only for values |fc| < k c where the cut-off k c is much 
smaller than the Fermi momentum k c <C kp. This im- 
plies that for not too strong interaction the ground state 
and the low energy excited states have negligible admix- 
tures of holes deep in the Fermi sea and particles with 
momenta \k\ — kp 3> k c . Therefore Tomonaga linearized 
the nonrelativistic quadratic dispersion in the vicinity 
of the two Fermi points ±kp. He then used Bloch's 
method of sound waves^ now called "bosonization" to 
obtain the ground state of the interacting system as well 
as its low energy excitations. As an application Tomon- 
aga calculated the density response of the (spinless) sys- 
tem in the long wavelength limit. This result was ob- 
tained before the pioneering paper by Bohm and Pines 3 
in which the authors introduced the "random phase ap- 
proximation" (RPA) for the (linear) density response of 
the (three-dimensional) interacting electron gas in order 
to obtain e.g. the collective plasmon mode. The RPA 
is a good approximation in the high-density limit. Var- 
ious alternative derivations of this approxmation can be 
given^, e.g. using diagramatic techniques or linearizing 
the time dependent Hartree-Fock (TDHF) equation^ and 
neglecting the exchange contribution and replacing the 
HF-eigenvalues by the noninteracting ones. 

If these additional approximations in the TDHF equa- 
tion are not made the corresponding response functions 
are called RPAE, where the "E" stands for exchange. 
While the RPA result is a simple expression in terms of 
the response function of the noninteracting system, the 
calculation of the RPAE response function usually re- 



to use TDHF which is a conserving approximation^ as 
the starting point to calculate the damping of the col- 
lective mode(s). We show that it is possible to solve 
the RPAE equation for the Tomonaga model analyti- 
cally when the two-body interaction v(q) takes a con- 
stant value for \q\ < k c and is zero otherwise. This exact 
solution reduces to the RPA result with the noninteract- 
ing response function ( "bubble" ) for strictly linear energy 
dispersion but differs markedly for nonlinear dispersion. 
The continua are in different places and the dispersion of 
the collective modes differ. This is most drastically seen 
for the collective spin mode which exists in RPAE but 
not in the bare RPA. 



II. LINERARIZED TIME DEPENDENT 
HARTREE-FOCK APPROXIMATION 



We want to describe the density response of a system of 
fermions on a ring of circumference L. The Hamiltonian 
is given by 
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fe 4 C L C L C fe4 C fc3 



(1) 

The index fc in addition to the momentum (possibly) in- 
cludes a spin index. The first term presents the kinetic 
energy with the nonrelativistic dispersion e k = k 2 /2m. 
The two-body matrix element v kl>k2>k3lki contains a fac- 
tor 8k 1 +k 2 ,k 3 +ki due to total momentum conservation. 
The remaining fc-dependences are specified later. 

In order to study the linear charge or spin density re- 
sponse the Heisenberg equations of motion for (c k c k ')t 
have to be solved in the presence of a time dependent 
external one-body potential with matrix elements Vkk' ■ 
Using the Hartree-Fock approximation to factorize the 
four fermion expectation values and linearizing the ex- 
pectation values around the HF ground state 



(clc k ,) t =6 kk ,f(e% F )+5(clc k ,) t 



yields 

d 
dt 



(4 c fe')t 



(2) 



(3) 
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/(efe' F )) Wk'k + v fc'fc4fcfc2^(ci 4 c fc2>t] . 



k 2 ,k 4 



where the HF-energies are given by 



(4) 



and v kltk2>kstki = v kl ,k 2 ,k 3 ,ki ~ v kl ,k 2 MM are the anti- 
symmetrized two-body matrix elements. At T = the 
Fermi function is a step function not depending on the 
interaction f(e k F ) = Q(k F - |fc|). Eqs. © are the lin- 
earized TDHF equations for the expectation values^ 

As we are interested in the long wavelength response 
we put k — > k a and k' — > k a + q, where \q\ <C kp and k a 
with a — R, L is in the neighborhood of the right or left 
Fermi point ±kp. Fourier transforming Eq. ([3]) yields 
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Vk, 



+q,k a 



■ ^ V k a +q,&,k a ,k'+q 
W 



,5(cl,C k > +g ) u 



where uj = to + iO. With the assumption \q\ kp the k' 
sum in Eq. ([5]) only can take values in the neighbourhood 
of the two Fermi points. This implies that two different 
types of two-body matrix elements contribute. For k' in 
the neighborhood of the same Fermi point as k a 

Vk a +q,a;k' a a';k a: cr;k' a +q,cr' = ~ [fif4(?) ~ 5<Ttr'94,(k a ~ K)] • 

(6) 

Here we have used the usual "g-ology" nomenclature^ 
and explicitely introduced the spin dependence. For the 



two-body matrix element with k a and k' w fc s where 
a differs from a no exchange contribution occurs using 
Tomonaga's assumptions as this would involve a momen- 
tum transfer of order 2kp, for which the two-body inter- 
action is assumed to vanish. Thus 

Vk a +q,cr;k' d a';k a ,cr;k' d +q,cr' = jT Q^io)- CO 

We separate the two contributions in Eq. (|5|) and ex- 
press the matrix elements of the external potential as 

V ka+q ^. ka ^ = V a , a {q)/L 



S ( c L,<T°ka+q,<r) 



1 /« F ) - f(ef! +q ) 



(8) 



+ ^4(9) - S a ^g 4 (k a - k' a ))S{cl latiT ,c k > a+qt(T >) u ] . 

k'.rr' 

For 32(9) = 34(9) = v(q) this are the RPAE equations for 
the Galilei invariant Tomonaga model. Because of the 
exchange term in the (^-interaction they usually have to 
be solved numerically^ They simplify considerably for 
the special choice for the g m (q), (m = 2,4) mentioned in 
the introduction 

9m(q) = g m Q{k c - \q\) ■ (9) 
With this assumption only expectation values of the type 

5{Pq,a,°)w = Yl S ( C k> a ,<T C k' a +<l><>)" ( 10 ) 

K 

and k' a — > k' a occur on the rhs of Eq. |(5J| . In Luttinger liq- 
uid terminology this are expectation values of the Fourier 
components of the densities of the right or left movers. 
Therefore Eq. {8} , after summing over k a , yields a closed 
set of equations for the 8{p q . a ^). Dropping the index u 
it reads 

S(Pq.a.a) = #H F W ) \V",o. (?) + 9i(^{Pq,a) ~ 8 {pq,a,<r)) 
+ 92S(pq,a)], (11) 



with 



«hf(«,w) 



1 V- MP - f($I + q) 



HF _ ,HF\ 
k a +q k a > 



(12) 



and the 8{pq^ a ) are the charge densities for the right and 
left movers 



With 



R a (q,u) 



1+94Rhf(i^) 



(13) 



(14) 
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Eq. dTTJ> reads 

S(Pq^a) = R a (q,u) V a ,a + giS{p Qja ) + g25{pq,a) (15) 

This equation is the starting point for the calculation of 
the charge and spin response. 

The equations for the charge response are obtained by 
summing Eq. (|15[) over the spin index 

(1 - 2g 4 R a )6( Pq , a ) - 2~g 2 R a 5{ Pq - a ) = 2R a V^ (16) 



with V4 (c) = J2a V<T, a /2. The solution is given by 



2R a (l - 2.g 4 i? a )K (c) + Ag 2 R a R a V^ 
(1 - 2g 4 R a )(l - 2g 4 R a ) - AglR a R« 
= R aa (q,uj)V^ +R a - a (q,uj)V^ c) (17) 



K<0 



The spin response 8{a q ^ a ) = J2a a &{pq,a,o) is obtained 
by taking the difference in Eq. (fT5|) . With Va^ = 
J2 a < J Vo-,a/2 one obtains 



S(a q , a ) u = 2R a = 



2i% F (g,^) 

1 +.94-RhfO?' w ) 



. (18) 



For the spinless model originally studied by Tomonaga 
the two gi terms on the rhs of Eq. (jTTj) cancel and one 
obtains^! 

si \ _ ffi-bja + .92-Rhf^hf^5 / 1Q \ 

°\Pq,a/u — , -2na na ' V 1J 7 

1 .y2 rl HF- n 'HF 

Equations (flT]) . (IS]), and (fT!)]) are the central results of 
this paper. The general behaviour and limiting cases of 
the corresponding response functions are discussed in the 
following section. 



the interaction and and therefore only differences of HF- 
energies occur in Eq. (Ti"2"|) . 

For \q\ < k c the behaviour of the response functions 
can be expressed in terms of two types of dimensionless 
quantities using m = k F /v F . 
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2ttv f 
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q 

2k F 



(21) 



If we define s a = 1 for a = R and s a — — 1 for a = L the 
HF response functions in Eq. (fT2"|) are given by 



Rn F (q,uj) 



1 1 

2~P2wv F 



log 



1 + «4 + (3 - s a Co 



1 + Q!4 — j3 — S a UJ 



> (22) 



with Co = uj / (v F q). 

The linear dispersion limit corresponds to the high den- 
sity limit /3 — > (fci? — > oo) 



/3 -> 



^hf(9,^) 



1 



2ttv f u> — s a (l + a^) 



(23) 



In this limit the functions i? a reduce to the noninteract- 
ing ones for linear dispersion 



R a (q,to) 



q 



2n cu — s a v F q 



(24) 



The well known results for the charge and spin mode 
of the Tomonaga model for linear dispersion 8 - ' 18 ' then 
follow from Eqs. (JTTj) and £jjl| as 



,2 

g,c 
2 

q,s 



= (v F qf [(1 + 2a 4 ) 2 - (2a 2 ) 2 ] 

= (v F q) 2 . (25) 



For the spinless model Eq. ([19]) yields the same form for 
the charge mode as in Eq. (|25|) but with 2ai — > a%. 



III. GENERAL RESULTS AND LIMITING 

CASES 

In this section we discuss the RPAE charge and spin 
response restricting ourselves to zero temperature. The 
basic building blocks of the RPAE response functions are 
the functions R a (q,uj) which via Eq. (fT4"]) are expressed 
in terms of the R^ SF (q,u). In order to calculate these 
functions we first have to determine the HF-eigenvalues 
e^? F defined in Eq. ((4]). For the special choice of the 
two-body matrix elements in Eqs. ©, ([7]), an d © one 
obtains for k = ±k F + k with |fc| < k c and k c < k F /2 



HF 

±k B 



± — k + const. 
2tt 



(±k F + k) 2 
2m 



±{v F + ^-)k+ — 
2tt 2m 



(20) 



c = ±v^ F k + — 
F 2m 



The constant is given by the Hartree contribution and 
the constant part of the Fock term. It is not needed 
here, as the Fermi function at T = is independent of 



A. The <?4 model 

While the spin response function in Eq. (|T5]) only de- 
pends on the 34 interaction, the charge response functions 
R aa i in Eq. (|17p involve a coupling of both Fermi points 
via the g 2 interaction . The chiral model obtained by 
putting g 2 = is called the 34 mode l 18 ' 19 . We discuss its 
response functions for the right movers and q > 0. The 
charge response follows from Eq. p7|) as 



08 -0) 



2R*(q,u) T - (c) 



l-2g 4 R R (q,u) 

2R§(q,Lu) Vl n 



Vj, c \q,u;). (26) 



l-2g i R^(q,to)i 



In the second line we also presented the result for the 
linear dispersion case. It looks tempting to generalize 
it to nonlinear dispersion by simply dropping the index 
"lin" and use instead the noninteracting response func- 
tion for quadratic dispersion which is obtained by putting 



a4 = in Eq. (122)) . As pointed out in the introduction 
this would be a bad approximation as it e.g. has the non- 
vanishing imaginary part of the response function in the 
wrong frequency range. To show this we use Eq. (TT4")) to 
express R R in terms of R^p 



2flg F (g,u;) 

1 -g 4 R^ F (q,u 



(27) 



For q > the continuous part of the imaginary part of 
this RPAE response function corresponding to the HF 
lp-lh excitations is different from zero for vpq(l + CV4 — 
P) < oj < vpq(\ + a>4 + (3), while the imaginary part 
of R R (q,ui) is different from zero for VFq{^ — P) < us < 
v F q{l+p). 

The RPAE charge mode follows from Eq. (gTJ) by 
putting the denominator to zero. The solution of the 
resulting transcendental equation is given by 



US 



RPAE 



VFq ( 1 + «4 + P coth — 

Ck4 



(28) 




FIG. 1: Dispersion ^^/(fF^F) of the collective charge and 
spin modes for the (74-model including spin for a repulsive in- 
teraction of dimensionless strength 04 = g4/(2nvp) = 0.5. 
The charge mode (full line) lies above the HF lp-lh contin- 
uum and the spin mode (dashed line) below it. This contin- 
uum is indicated by the hatched area. 



which in the limit P — > agrees with the 0:2 = result 
in Eq. (|25| . The mode is always above the the HF lp- 
lh continuum as also shown in Fig. 1. This is not the 
case for the RPA-charge mode which is obtained from 
1 - 2g A R R = as 



Lo^ A ^v F q[l + Pcoth^- 

2«4 



-q,c 



(29) 



It follows from Eq. (TTg)) that the RPAE collective spin 
mode is generally independent of §2 ■ Its location is given 
by 



US 



RPAE 
q,s 



VFq I 1 + «4 — /3coth 



fl'4 



(30) 



which in the limit P — > agrees with the result in Eq. 
(|25p which is independent of the interaction strength for 
the Tomonaga model with linear dispersion. The RPA 
spin response function 2R^(q,uj) has only a continuous 
non- vanishing imaginary part. This is different for the 
RPAE result for nonlinear dispersion, where the delta 
spike of the spin mode is always below the HF lp-lh 
continuum. 

In Fig. 1 we show the dispersion of the modes relative 
to the HF-continuum for a 4 = 0.5. 

The imaginary part of the RPAE spin response func- 
tion as a function of us for fixed P = 0.1 is shown in Fig. 
2, where it is compared to the result of the (bare) RPA 
result, which shows no collective mode. The weight of 
the RPAE collective spin mode which is slightly below 
u> = 1 carries 94.8 per cent of the total RPAE spectral 
weight. 



pi 
B 



CO/(v q) 



1.5 



FIG. 2: Imaginary part of the spin response function for the 
gi model for q/kF = 0.2 as a function of us for 04 = 0.5, 
as in Fig. 1. The full line shows the RPAE result, with the 
collective spin mode indicated by the arrow. The dotted curve 
shows the RPA result which only has a continous part but no 
collective spin mode. 



B. Galilei invariant model: 



32 = 34 



wo 



In this subsection we discuss the model that describes 
interacting fermions with a spin independent two-body 
interaction v(\xi — Xj\) on a ring with 



v(q) = v O(k c 



(31) 



and k c < kp/2. 

Then the expression for the total charge response in 
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0.1 



q/(2k F ) 
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FIG. 3: Dispersion of the normalized charge velocity 
v c (q)/vF = ui q , c /(qvF) for the Galilei invariant model (32=34) 
including spin for a repulsive interaction of dimensionless 
strength a = v /(2ttv f ) = 0.25. The RPAE result (full line) 
lies above 1 + a + f3 indicated by the dashed line, marking the 
top of the HF lp-lh continuum. The RPA charge velocity 
(dashed-dotted line) crosses the dashed line and only stays 
above 1 + f3 corresponding to the bare lp-lh continuum (not 
shown) . 



Eq. (fTT|) simplifies to 



R{q,u) 



09 — 0) 



1 - v R(q,u) 

Ro(q,u>)lin 
1 - V R (q,uj)un 



(32) 



where in R = 2(R R +R L ) and R = 2(R$+R£) the factor 
two is the due to the spin degeneracy. For the spinless 
model one obtains the same result but this factor of two 
is missing. 

As in the previous subsection it would be a bad ap- 
proximation to generalize the bare RPA in the second 
line in Eq. (f3"2"| (which provides the exact result for lin- 
ear dispersion) by simply dropping the index "lin" . 

The transcendental equation which results from 
putting the denominator in Eq. ()32p equal to zero can 
no longer be solved analytically. In Fig. 3 we show nu- 
merical results for v c (q) = w 9;C /|g| for the spinfull model 
for a = v q /(2ttv p ) = 0.25. 

As a comparison we also show the corresponding bare 
RPA result. For < a < 2 the charge velocity for lin- 
ear dispersion v c — vf^/I + 4a is larger than v^ F = 
v F (l + a). For these a values the RPAE charge mode 
is above the HF lp-lh continuum. For larger repulsive 
interactions a > 2 the charge velocity for linear disper- 
sion is smaller than v^ F which puts the RPAE charge 



mode below the HF lp-lh continuum. 

In the spinless model the RPAE charge mode for all 
a > lies below the HF lp-lh continuum. This can easily 
be seen by examining the denominator in Eq. (|19|) . This 
is qualitatively different from the RPA result which has 
the charge mode above the noninteracting continuum. 

In an approximation beyond RPAE the collective 
modes discussed above are expected to aquire a finite 
lifetime leading to a broadening of the delta peaks in the 
imaginary part of the response function. Because the 
higher HF p-h continua have an additional extent for en- 
ergies above the lp-lh continuum but not below it, the 
damping of the collective modes will depend crucially on 
how they are located relative to the HF lp-lh continuum. 



IV. SUMMARY 

In this paper we have revisited the old problem of the 
long wavelength linear charge and spin response of in- 
teracting fermions in one spatial dimension. As there is 
renewed interest in the questions of the lifetime of the 
ideal bosonic modes of the Tomonaga model due to de- 
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we 



viations from the linear electronic dispersion 
hope to have clarified the often quoted statement that 
"RPA becomes exact in the long wavelength limit" . We 
have shown that a naive generalizaton of this result fails 
badly. 

We have not presented a calculation for the lifetimes 
but only set the proper stage from where to start a per- 
turbative approach. In order to avoid tadpole diagrams 
HF-propagators have to be used. The use of bare prop- 
agators is useful only for the case of linear dispersion 
because of the "closed loop theorem" 

For the spinfull model and not too strong repulsive 
interactions the collective charge and spin modes are 
above respectively below the HF lp-lh continuum. As 
the higher particle-hole continua for fixed q with particle 
hole pairs around both Fermi points are higher in energy 
the damping of the spin mode is expected to be quite dif- 
ferent from the charge mode. For the spinless model and 
repulsive interaction the RPAE charge mode is always 
below the lp— lh continuum. Therefore we also expect 
the spinfull and the spinless model to behave differently 
concerning the damping of the charge mode for weak re- 
pulsive interactions. 
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